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Topological charged black holes coupled with a cosmological constant in R 2 x X D ~ 2 spacetimes 
are studied, where X ~ is an Einstein space of the form (D_2) i?As = k(D - 3)h A B- The global 
structure for the four-dimensional spacetimes with k — is investigated systematically. The most 
general solutions that represent a Type II fluid in such a high dimensional spacetime are found, and 
showed that topological charged black holes can be formed from the gravitational collapse of such a 
fluid. When the spacetime is (asymptotically) self-similar, the collapse always forms black holes for 
k — 0, — 1, in contrast to the case k = 1, where it can form either balck holes or naked singularities. 

I. INTRODUCTION 

Lately black holes in high dimensions have attracted a great deal of attention in the gravity-gauge theory correspon- 
dence [1], and been further promoted by theories of TeV gravity, in which high dimensional black holes are predicted 
to be produced in the next generation of colliders [2] . 

In black hole physics, one of the fundamental features is the topology of a black hole. In four-dimensional asymptot- 
ically flat stationary spacetimes, Hawking first showed that a black hole has necessarily a S 2 topology, provided that 
the dominant energy condition holds [3]. Later, it was realized that Hawking's theorem can be improved in various 
aspects, see [4] and references therein. However, once the energy condition is relaxed, a black hole can have quite 
different topologies. Such examples can ocurr even in 3 + 1 dimensional spacetimes where the cosmological constant 
is negative [5]. In high dimensional spacetimes, it was found that even if the energy conditions hold, the topology is 
still not unique. In particular, a five-dimensional rotating vacuum black hole can have either a S 3 topology [6], or a 
S 2 x S 1 topology [7], although a static charged dilaton black hole in high dimensional asymptotically flat spacetimes 
indeed has an unique topology, which is a (D — 2)-sphere [8]. When the cosmological constant is different from zero, 
similar to the four-dimensional case, static black holes can have different topologies [9]. 

In this paper we shall study the formation of topological black holes from gravitational collapse of a type II fluid 
in high dimensional spacetimes. Specifically, the paper is organized as follows: In Sec. II we present a general 
R 2 x X D ~ 2 decomposition. In Sec. Ill, assuming that X D ~ 2 is an Einstein space with a constant curvature, 
(-°~ 2 ) Rab = krjtiAB, we re-derive the charged solutions coupled with a cosmological constant in any dimensional 
spacetimes, without assuming that the spacetime is static. In Sec. IV we systematically study the global structure 
for the case D = 4 and fcrj = 0, while in Sec. V, all the type 77 fluid solutions in D-dimcnsional spacetimes are given. 
In Sec. VI we study the formation of topological charged and uncharged black holes from the gravitational collapse 
of such a fluid, while in Sec. VII our main conclusions are summarized. There is also an appendix, in which trapped 
surfaces and apparent horizons are defined. 

Before proceeding further, we would like to note that the formation of topological black holes from gravitational 
collapse in four-dimensional spacetimes was studied in [10], while gravitational collapse in high dimensional spherically 
symmetric spacetimes was investigated in [11]. 
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II. R 2 x X u - 2 DECOMPOSITION 

In this paper let us consider a D-dimensional spacetime described by the metric 

ds 2 = g^dx^dx" = j ab (x c )dx a dx b + s 2 (x c )h AB (x c )dx A dx B , (2.1) 

where we use Greek indices, such as, /it, v, A, to run from to D — 1, lowercase Latin indices, such as, a, b, c, to 
run from to 1, and uppercase Latin indices, such as, A, B, C..., to run from 2 to D — 1. Clearly, the above metric 
is invariant under the coordinate transformations, 

x a = x a (x /b ), x A =x A (x' B ). (2.2) 

Introducing the quantities -f ab (x c ) and h AB (x c ) via the relations, 

7 ac 7cfc = SI h Ac h CB = S A , (2.3) 

we find that 

Then, the non-vanishing Christoffcl symbols are given by 

(D)-pa _ (2) -pa (D)pa _ _^,„ S o- a 

where s j0 = ds/dx a , s' a = j ab s, b , ^T^, and are the Christoffel symbols calculated, respectively, from j ab (x c ) 

and h (x ), and n = D — 2. From Eq.(2.5) we find that the Ricci tensor, defined by, 

has the following non-vanishing components, 

^R ab = ^R ab -—V a V b s, 
s 

MRab = {n) R AB - [sUs + (D - 3) (V c s) (V c s)] ^b, (2.7) 

where □ = 7 ah V a V{,, and V a denotes the covariant derivative with respect to j ab . 

To calculate ^R ab , we first note that the Riemann tensor ^R a bcd has only one independent component, say, 
( 2 )i?oioi- Then, in terms of ^ 2 ^i?oioi we have 

(2) R abcd = (2) iW {Ol8% - 6°jl6l6° d - Sl5 b S A + SfflSlft) , 
^R ab ee ^R\ cb = ^R 0101 (^5y b - 7 01 (6°Jl + Sy b ) + 7 00 <W • (2.8) 

On the other hand, for the second block of the metric, we shall consider only the case where \iab represents a 
n-dimensional surface with constant curvature, 



(n) 



Rab = k D h AB (x c ), (2.9) 



where, without loss of generality, one can always set k B = (D — 3)k with k = —1, 0, 1. Then the surfaces of constant 
x a will be referred, respectively, to as hyperbolic, flat and elliptic. 

III. TOPOLOGICAL CHARGED BLACK HOLES COUPLED WITH THE COSMOLOGICAL 

CONSTANT 

Using the gauge freedom of Eq.(2.2), in this section we shall choose the coordinates such that 

ltr (x c ) = 0, s(x c ) = r, (3.1) 
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where x c = {t, r}. This will be referred to as the Schwarzschild gauge. Then, we find that 

ds 2 = ^{t.r) dt 2 + e 2*(t,r) dr 2 + ^ h^C^A^B 

For such a form, there still remains the gauge freedom 

t = t(t). 

Later, we shall use it to fix some integration functions. For the metric (3.2) we find 

= 2 {e- 2 * [* „ + * t (* t - $ t )] - e- 2 * [* rr + * r ($ r - * r )]} . 

It can be also shown that 

Va V fc r = -e^^,Sy b *, t (Sy b + 6ffl) 
Ur = e- 2 ^ (* r -* r ). 

Substituting Eqs.(2.9), (3.4) and (3.5) into Eq.(2.7) we obtain 

MRoo = - [*, a + *, t (*,t-*,t)] 



(3-2) 
(3.3) 



(3.4) 



(3.5) 



2(*-*) 



+e 



$ rr + $ r $ r - * r + 



L> - 2 



Wi2 n = e 2 (*-*) [*,« + *, t (*,*-*,*)] 

/ D - 2 



= {k D - [r ($ r - tt, r ) + (D - 3)] e- 2 *} /i AB . 
On the other hand, the D-dimcnsional Einstein-Maxwell equations with the cosmological constant, Ad, read 

{D) R^ - \g»» {D) R = 8ttG d We^ - A D9 ^, 

where the energy- momentum tensor ^E^ is given by 

(D) E - 1 ( (D) p (D) p ° _ 1 (D)p (D)p a f)\ 

~ AttG d V 4 / 

with the electromagnetic field ^F^ satisfying the Maxwell equations, 

D v (D) F^=0, 



(3.6) 
(3.7) 

(3.8) 
(3.9) 

(3.10) 
(3.11) 
(3.12) 



where D v denotes the covariant derivative with respect to the D-dimensional metric <7 M „. From the symmetry of the 
spacetime, without loss of generality, we can always assume that ^F^ v is function of t and r only, and takes the 
form, 



Substituting Eqs.(3.11) and (3.12) into Eq.(3.10) we find that 



(3.13) 



^R^=F 2 2 



+ 



2A D 
D-2 



D-2 



(3.14) 



Combining Eq.(3.7) with the above equation, we obtain 

#(t,r) = tf(r). (3.15) 
On the other hand, from the Maxwell field equation (3.12) we have 

F(t,r) = ^e^, (3.16) 

where Co is an integration constant. Inserting Eqs.(3.15) and (3.16) into Eqs.(3.6)-(3.9), we find the following 
independent equations, 

D-2. T . d 2A D 



<V + *,r (*,r - *,r) " — r = ^ ^) ■ (3-18) 

*„e» - r (*, r - *,) - (D - 3) = e» ( (D _ 3 ^ (D - 3) + ^r>) , (3-19) 

where Ci = 2(D - 3)C$/(D - 2). From Eqs.(3.17) and (3.18) we obtain 

*(t,r) = -*(r) + $o(t), (3-20) 

where $o(i) is an arbitrary function. However, using the coordinate transformation (3.3), we can always set $o(i) = 0, 
a condition that will be assumed below. Then, substituting the above expression into Eq.(3.19) and integrating it, we 
obtain 

* = -* = H fc -^ + ^-A- 2 ), (3 ' 21) 

where 

A5 (D-^-2V (3 - 22) 

and M and Q are two arbitrary constants, related to the total mass Mu and charge Qd of the spacetime via the 
relations [12], 

where As denotes the total area of the (D — 2)-dimensional surface spanned by Hab, 

A s = J Vhd°- 2 x. (3.24) 
It can be shown that the solution (3.21) also satisfies Eqs.(3.17) and (3.18), with the electromagnetic field given by 

(D) ^ = { {D ~ 2 f - 3) ) ^ (Wl - ^) . (3.25) 

When the surfaces of constant t and r are not compact, clearly A s becomes infinity, so do M D and Qn- However, 
one can still keep the quantities m and q finite, where 

M D Qd OR , 

m=— , 9=^-, (3.26) 

which may be interpreted as the mass and charge per unit area. 

On the other hand, from the above derivation of the solutions one can see that they represent the most general 
solutions of the Einstein-Maxwell field equations coupled with the cosmological constant in D-dimensional spacetimes 
described by the metric (3.2). In general the spacetime is singular at r = and all the other singularities are coordinate 
ones. This can be seen, for example, from the Ricci scalar, which now is given by 

<"'«=- 4)02 +^A»- 

Thus, to have a geodesically complete spacetime one needs to extend the spacetime beyond these coordinate singu- 
larities. In the next section, we shall restrict ourselves only to the case where kn = and D = 4. 
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IV. FLAT TOPOLOGICAL BLACK HOLES IN FOUR-DIMENSIONAL SPACETIMES 



To study the global structure of the solutions given in the above section, in this section we shall restrict ourselves 
only to the flat case, ko = 0, and assume that the spacetime has only four-dimensions, D = 4. Then, the surface of 
constant t and r can be a plane, a cylinder, a Mobius band, a torus, or a Klein bottle, depending on how to identify 
the two coordinates in X 2 [13]. 

Setting k = and D = 4 in Eq.(3.21) we find that the corresponding metric can be written in the form, 

ds 2 = -f(x)dt 2 + f-\x)dx 2 + x 2 (dy 2 + dz 2 ), (4.1) 

where 



... 26 q 2 A 2 
f(x) = — + — - -x\ 
x x z 3 



with b = —M, q = Q and A actually being A^. To consider the most general case, we assume that the mass parameter 
b and the coordinate x take their values from the range — oo < b, x < +oo. It is interesting to note that if we make the 
replacement (M, r) — > (— M, —r), we find that the Reissner-Nordstrom solution coupled with a cosmological constant 
remains the same. That is, the physics for M > 0, r > is the same as that for M < 0, r < 0. When the spacetime 
has plane symmetry, the range of the coordinate x is — oo < x < +oo. This is the main reason why here we are 
allowed ourselves also to consider the case where b > 0. 

It should be noted that the above solutions was also studied [5]. Since here we would like to give a systematical 
study for this case, some of the materials to be presented below will be unavoidably overlapped with some presented 
there. 

From Eq.(4.1) we can see that the metric is singular on the hypersurfaces f(x) = 0. However, these singularities 
are coordinate singularities, except for the one located on the hypersurfacc x = 0. This can be seen, for example, 
from Eq.(3.27) as well as the corresponding Krctschmann scalar, which now is given by 



6(q 2 + bx) 2 + q i ]+^A 2 \. (4.3) 



The above expression shows clearly that the Kretschmann scalar becomes unbounded only when x — > 0. As \x\ — > +oo 
the spacetime is asymptotically dc Sitter or anti-de Sitter, depending on the signs of the cosmological constant, A. 
Thus, only the singularity at x — is a spacetime curvature singularity, which divides the whole (t, a;)-plane into 
two causally disconnected regions, x > and x < 0. All the other singularities are coordinate ones, and to have a 
geodcsically complete spacetime, we need to extend the metric beyond these singularities. 

In the stationary axisymmetric case, a significant quantity that characterizes a horizon is the surface gravity [14]. 
In the following we shall derive this quantity for the spacetimes described by the metric (4.1) [15]. To this end, let us 
consider the four-acceleration of a static observer with the four- velocity u M = / _1 / 2 ^, 

a" = uP; v u v = k(x)6£, (4.4) 

where 

K(x) = ^f, x {x). (4.5) 

From the above equations we can see that in order to stay at a fixed point x, the observer has to fire a jet engine whose 
thrust per unit mass is |«(a;)| in the positive x— direction for f, x (x) > 0, and negative x— direction for f, x (x) < 0. 
On the horizon x = x , where f(x ) = 0, the quantity k defined by 

k = k(x ) = (x ) , (4.6) 

is called the surface gravity of the horizon 1 . 



J Note the similarity to the spherical case: ds 2 = f(r)dt 2 — f(r)~ 1 dr 2 — r 2 (d9 2 + sin 2 6d<f> 2 ). For a static observer with the 
four-velocity u M = /~ 1//2 <5f, its four-acceleration is given by a M = = (f, r /2)5p. The quantity k = f, r (r g )/2 is called 

the surface gravity of the horizon at r — r g [14]. 
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To extend the spacetime across the horizon, we can first introduce a new coordinate x* by [16] 



x = 



dx 



/(*)' 



and then two null coordinates, v and w, via the relations 



v = e K(t+x *\ w = -e-^ t - x "\ 



In terms of v and w, metric (4.1) becomes 

ds 2 = -n- 2 f(x)e- 2KX *dvdw + x 2 (dy 2 + dz 2 ). 
To study the global structure of the spacetime, it is found convenient to distinguish the following cases: 



A) 0, q ^ 0, A = 
C) b= 0, q ^ 0, A ^ 
E) 0, q = 0, A = 
G) 6= 0,q = 0, A ^ 0. 



B)b^0, q = 0, A t^O; 
^)6^0, g ^0,A^0; 
F) 6 = 0, q ^ 0, A = 0; 



(4.7) 

(4.8) 
(4.9) 



(4-10) 



In the above classification, Case E) corresponds to the Taub vacuum solution [17] and the properties of which are 
well-known. In particular, it has a naked singularity at x — 0. Case F) represents a pure electromagnetic field with a 
naked singularity at x = 0, while Case G) corresponds to the de Sitter or anti-de Sitter spacetime, depending on the 
sign of the cosmological constant. The global structure of the latter is studied extensively in [18]. Therefore, in the 
following we shall focus our attention only on the first four cases. 



A. 6^0,g^0,A = 0. 

When A = 0, we find 




(4.11) 



FIG. 1. The function f(x) given by Eq.(4.11) in the text, (a) The case where b < 0, in which the function f(x) is zero at 
x — xo = q 2 /(2\b\) and has a minimum at x = xi = 2xq, where /i(a;i) = —b 2 /q 2 . (b) The function f(x) for b > and now wc 
have xo = -q 2 /( 2 &) < and fi(2x ) = -b 2 /q 2 . 



G 



Case A.l) b < 0: In this case the funtion f(x) is positive in the whole region x < [cf. Fig. 1(a)], and the 
singularity at x — is naked. The corresponding Penrose diagram is that of Fig. 2. Note that now k{x) = f'(x)/2 = 
(\bx\ + g 2 )|a;| -3 is also positive in this region, so the naked singularity at x = produces a repulsive force to a 
stationary observer who stays in this region. Consequently, the naked singularity should effectively have negative 
gravitational mass. 




FIG. 2. The Penrose diagram for the solution of Eq.(4.11) with b < in the region x < 0. The singularity at x = is timelike 
and naked. The vertical curves are the hypersurfaces x = Const, and the horizontal dotted lines are the ones t = Const. 



However, in the region x > 0, the function f(x) becomes zero at the point 



x = |^>0, (6<0), (4.12) 
and the corresponding surface gravity is given by 

« = -^, (6<0). (4.13) 

Thus, a stationary observer near the hypersurface x = xq in this region feels a repulsive force from the region 
< x < xo, and the totally effective gravitational mass should be also negative in this region. Integrating Eq.(4.7) 
we obtain 



On the other hand, the coefficient g vw now is given by 



g vw = -(V2kx)- 2 cxp (-^ (q 2 + \b\xf) , (4.15) 

which shows that the metric now is well defined across the hypersurface x = x . Combining Eqs.(4.8) and (4.14), we 
also have 

vw = - (q 2 - 2\b\x) exp {^-x 2 + Hi*U) , (4.16) 
v= (q 2 -2\b\x) 1/2 exp(-^t+^x 2 + ^x) (Region!). (4.17) 
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From the above expressions we can see that the coordinate transformations (4.8) are restricted only to the region 
where v > and w < 0, which will be referred to as region / [cf. Fig. 3]. To extend the metric (4.9) to cover the 
whole spacetime, one simply takes the range of v and w as — oo < v, w < +oo. Then, in the (v, w)— coordinates we 
obtain three extended regions /', II and II', which are absent in the (t, x)— coordinates. In each of the three extended 
regions, the coordinate transformations from (v,w) to (t,x) are given by 

v = -(q 2 - 2\b\x) 1/2 cxp + ^ + , (Region 1% 

v = (2\b\x - q 2 ) 1/2 cxp (-^-t + ^x 2 + ^x] , (Region II), 

\ q 4 <r <r ) 

v = -(2\b\x-q 2 ) 1/2 cxp(^-^t + ^x 2 + ^x^ , (Region II'), (4.18) 

where w can be found through Eq.(4.16). Then, the corresponding Penrose diagram is that of Fig. 3. From there we 
can see that the hypersurfaces x = xq (or vw = 0) are Cauchy horizons that separate the two regions / and I' from 
the two asymptotically flat (only in the x— direction) ones II and II' [18]. The time-like coordinate t is past-directed 
in region / and future-directed in region as we can see from Eqs.(4.16)-(4.18). Across the horizons, t becomes 
space-like, while x becomes time-like. 




FIG. 3. The Penrose diagram for the solution of Eq.(4.11) with b < in the region x > 0. The singularity at x = is timelike 
and the lines be and ad where x = xo are Cauchy horizons. 



Case A. 2) b > 0: In this case the spacetime properties can be obtained from the case b < by the replacement 
(b, x) — ► (—6, —x). In particular, in the region x > the function f(x) is strictly positive and the singularity at x = 
is naked. The corresponding Penrose diagram is similar to that given by Fig. 2, but now with x > 0. In the region 
x < 0, the metric becomes singular at the hypersurface x = xq = —q 2 /(2b) [cf. Fig. 1(b)], on which we now have 

46 3 

K = — > 0, (b > 0). (4.19) 

That is, a stationary observer near the hypersurface x — xq in this region feels a repulsive force from the region 
xq < x < 0. After extending the spacetime beyond this surface, the corresponding Penrose diagram is similar to that 
of Fig. 3 but with x < 0. 



B. 6^0,g = 0,A^0. 



When g = 0, we have 



, A f 6b 



(4.20) 
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In this case it is found convenient to distinguish the four subcases, 

1) A > 0, b > 0; 2) A > 0, b < 0; 

3)A<0, 6<0; 4)A<0, b>0. (4.21) 

Case B.l) A > 0, b > 0: In this subcase, we find that f(x) is negative for all x < 0. Thus, in this region x is 
timclikc and t is spacelike. Then, the singularity at x — is spacelike and naked. So, in this region the spacetime has 
a global structure that is quite similar to some cosmological models. 



f(x) 



Xi 


l\ 












fi 




X 



FIG. 4. The function f(x) given by Eq.(4.20) in the text, where Xl = -(3fe/A) 1/3 and = -[(3&) 2 A] 1/3 < 0. 



In the region x > the function f(x) is positive only for < x < x n [cf. Fig. 4], where 




(4.22) 



On the hypersurface x = xq we have f(xo) = 0, which, as shown before, is only a coordinate singularity. Thus, we 
need to extend the solution beyond this hypersurface. Before doing so, we first note that now we have 

k = -^x o <0, (4.23) 

that is, an observer near the horizon will feel a repulsive force from the region < x < x . Then, we would expect 
that the singularity at x = has an effectively negative gravitational mass. This is understandable, since in the 
present case we have A > 0, which is energetically equivalent to a matter field with a negative mass density. 

Following the same line as outlined in the last subcase, we find that the corresponding Penrose diagram is given by 
Fig. 5. The hypersurfaces x = x (or vw = 0) are Cauchy horizons. The singularities at x = now are also timelike. 
The coordinate t is past-directed in region / and future-directed in region 
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FIG. 5. The Penrose diagram for the solution of Eq.(4.20) with b > and A > in the region x > 0. The singularity at 
x — is timelike and the lines be and ad where x = xo are Cauchy horizons. On the curved lines ab and cd we have x — oo. 



Case B.2) A > 0, b < 0: From Eq.(4.20) we can see that this case can be obtained from the last case by the 
replacement (x, b) — > (—x, —6). Thus, the spacetime structure for this case can be also obtained from the ones given 
in the last case by this replacement. In particular, now the solution in the region x > may represent a cosmological 
model, and the spacetime has a naked singularity at x = 0, which is spacelikc. In the region x < 0, the corresponding 
Penrose diagram is that of Fig. 5 but now with x < 0. 




Case B.3) A < 0, b < 0: In this case, we find that 

/(*) = 
/'(*) = 



3x 
\2A\ 
3x 2 



66 

1 
36 
A 



(4.24) 



Thus, in the region x < the function is always positive as shown by Fig. 6. As a result, the singularity at x = is 
naked. On the other hand, for any hypersurface x = C, its normal vector n a is given by n a = d(x — C)/dx a = 5*. 
Thus, we have 



+oo, 



(4.25) 



as an — oo, which means that the spatial infinity x = — oo in the present case is timelike. Then, the corresponding 
Penrose diagram is given by Fig. 7. 

It is also interesting to note that now the acceleration of a static observer is given by Eq.(4.4) with 



«(*) = ^ ( ' + 



< 0, x < xi, 
= 0, x = xi, 
> 0, xi < x < 0, 



(4.26) 



where x\ 



-|36/A| 



1/3 
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b 

FIG. 7. The Penrose diagram for the solutions given by Eq.(4.24) in the region x < for A < and b < 0. 



In the region x > 0, we have 



fix) 



< 0, x < x , 
= 0, x = x , 
> 0, x > xo, 



where x = |6fo/A| 1 / 3 , and 



K = \f'{ x o) 



x > 0. 



(4.27) 



(4.28) 



Extending the solutions into the region < x < Xo, we find the following relevant quantities 

2(x 2 + x Q x + xlf/ 2 f r- {2x + x Q 

9vw = oi a i 2 exp <^ -V3 arctg 



vw = 



i\K\xlx 

X — Xq 



V = 



(x 2 + X X + Xq) 1 / 2 

(x-x a y/ 2 



exp 



V3 arctg 



(x 2 + x x + a;§)V4 
exp 



A V3 /2a; + a;o s 

2^ + - arctg l^7. 



•s/3x 
2a; + x 
V3xo 



(Region I). 



(4.29) 



(4.30) 



The corresponding Penrose diagram is given by Fig. 8, which shows that now the hypersurfaces x = x (or vw = 0) 
represent event horizons, and the singularities at x — are space-like. Regions iT and //' are two catastrophic regions, 
while regions / and /' are two asymptotically anti-de Sitter regions. The time-like coordinates t is future-directed in 
region I and past-directed in region Region II can be considered as a black hole, while region IV as a white hole. 
Since now the surface gravity of the horizon is positive, the spacetime singularity at x = is gravitationally attractive 
and expected to have an effectively positive gravitational mass. 
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x = 

FIG. 8. The Penrose diagram for the solution of Eq.(4.24) with b < and A < in the region x > 0. The singularity at 
x = is spacelike and the lines be and ad where x — xo represent event horizons, while the curved lines ac and bd represent 
the spatial infinity x — oo. 



Case B.4) A < 0, b > 0: This is the case which corresponds to Case B.3). As a matter of fact, replacing x by —x, 
we can obtain one from the other. Therefore, the corresponding Penrose diagrams in the two subcases are the same, 
after the above replacement is taken into account. 



C. 6 = 0,g^0,A^0. 

When 6 = 0, from Eq.(4.2) we can see that f(x) = f(—x). Thus, the spacetime has reflection symmetry with 
respect to the hypersurface x = 0. Then, without loss of generality, in the following we shall only consider the region 
x > in this case. 

When A < 0, f(x) is always greater than zero for any x. Consequently, the singularity at x = is naked, and the 
corresponding Penrose diagram is similar to that of Fig. 7. 
When A > 0, we have 

f(x)^^(xl+x^(xl-x"), (4.31) 

but now x a is given by 

1/4 

(4.32) 

From the above equations we can see that f(x) > holds only in the region < x < x . On the hypersurface x = x a 
the surface gravity is given by 

2A 

k = -—x <0, (4.33) 

that is, the effective gravitational mass in the region < x < xa should be negative. One can show that the 
corresponding Penrose diagram is similar to that of Fig. 5. 




D. &/0,g/0,A/0. 

When none of the three parameters are zero, we need to distinguish the four subcases as denned by Eq.(4.21). 
However, as shown in Cases A) and B), the properties of the spacetimes for the case b < can be obtained from those 
for the case b > by simply replacing x by —x. Thus, we actually need to consider only the subcases A > 0, b > 0, 
and A < 0, b < 0. 

Case D.l) A > 0, b > 0: In this case, it can be shown that f(x) takes the form 

f(x) = ^{x+ - x)(x - x_)(x 2 + Bx + C), (4.34) 
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where 



B = x++x_, C 



■ X+X- 



and x± are the two real roots of the equation f(x) — 0, given via the relations 

66 „ 3q 2 



B{x\ + x 2 _) = ^, Cx + x- 



A 



(4.35) 
(4.36) 



with the property 

x + > 0, x- < 0. (4.37) 

From the above equations we can see that f(x) > is true only in the region X- < x < or < x < x + . Then, we 
need to extend the solutions beyond both of the hypersurfaces x — x±. 

Let us first consider the extension across the hypersurface x = x + . We first note that 

n\ x=x+ = - g^-(z+ - x-){x* + + Bx + + C)<0. (4.38) 

Thus, a static observer near the horizon will feel a repulsive force from the region < x < x+, which means that the 
effective gravitational mass in this region is negative. One can also show that the Penrose diagram is similar to that 
given by Fig. 5. 

On the other hand, the surface gravity at x — x_ is given by 



^2"(x + - x-)(x 2 _ + Bx_ + C) > 0, 



(4.39) 



which means, similar to that in the region x > 0, now a static observer near the horizon will also feel a repulsive force 
from the region < x < 0. It can be also shown that the corresponding Penrose diagram is that of Fig. 5, but with 
x < 0. 

Case D.2) A < 0, b < 0: In this case we find that it is convenient to distinguish the following three subcases: 

i) \q\ < q c , ii) \q\ = q c , Hi) \q\ > q c , (4.40) 

where 



|A| 



-1/6 



2/3 



Case D.2.i) \q\ < q c : In this case, it is found that the function f(x) can be written as 

f{x) = ^(x - x+)(x - x_)(x 2 + Bx + C), 

where B, C and x± are given by Eqs.(4.35) and (4.36), but now with the properties 

x+ > 0, x_ > 0, (\q\ < q c ). 

From Eqs.(4.42) and (4.43) we can see that in the region x < the function f{x) is positive, and the singularity at 
x = is naked [cf. Fig. 9]. The corresponding Penrose diagram is similar to that given by Fig. 7. 

f(x) 



(4.41) 

(4.42) 
(4.43) 





f >\ 






\x_ 


/ x + - 





\ / X 
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FIG. 9. The function f(x) given by Eq.(4.42) in the text, for A < 0, b < and |g| < q c , where fi(xi) = [(36) 2 |A|] 1/3 > 0. 



When x > 0, f(x) is non-negative only in the regions x > x + or < x < x_. Thus, an extension across the 
hypcrsurfaces x = x± is needed. Let us first consider the extension across x = x + . 

The extension across the hypersurface x = x + : Following previous cases, we find the following, 



— ln(x - x+) + ^ ^{x - x+) n , 



K+ = 

m = 

9vw 



2k + 

f,x (X+) 

2 

x — x + 



n=l 



|A| 



— — (x + - x_)(x + + 6x + + c) > 0, 
ox i 



a n = 



1 cf l ff (x) 



g(x) ' n! dx™ 

-^(x)]" 1 exp \-2k + 2 — " *+)" 



!)io = — (x — x + ) exp 



u = (x — X+) 1 / 2 exp K + 



( QC 



n=l 



n 



-(x-x+r , 



n=l 



(Region I). 



(4.44) 



Then, we get four regions, 7, I', 77 and II' in the (v, w)-plane, as shown in Fig. 10. The hypersurfaces x = x + now are 
event horizons. The coordinate t is future-directed in region 7 and past-directed in region 7'. Note that as x — > +oo, 
we have f(x) — > +oo. Therefore, the hypersurface x = +oo is timclike and represented by the curves ac and bd in 
Fig. 10. 

The extension across the hypersurface x = x _: Now we find that 



1 °° b 

= ln(x — x-) + > — (x — x_) 

2k- ^ n 



m = 



2k 

f,x (x-) 
2 

x — X- 
S(x) 



n=l 



-7r^r(x+ - x_)(x 2 _ + 6x_ + c) < 0, 
6x1 



b n = 



1 d n g(x) 



n\ dx n 



g vw = -^nlgix)]- 1 exp \ -2k. £ ^(x - x_)" 



ra=l 



VW 



(x — x_) exp < — 2k_ ^ — (x — X-) n 



n=l 



(x — x_) 1//2 exp < k 



t+^(x- x _r 



n=l 



(Region III). 



(4.45) 



Clearly, by this extension, we get two more regions 777 and 777', as shown in Fig. 10. The singularities at x = are 
time-like and repulsive, since now we have K- < 0. The hypersurfaces x = x_ represent Cauchy horizons. Extensions 
across these hypersurfaces arc further needed in order to have a geodesically complete spacetime. However, since 
these surfaces are Cauchy horizons, the extensions across them are not unique, a situation quite similar to that of the 
Reissner- Nordstrom (RN) or Kerr solution [18]. 
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FIG. 10. The Penrose diagram for the solution of Eq.(4.42) with \q\ < q c in the region x > 0. The lines ad and be where 
x = x+ represent event horizons, and the ones gb, ha, cf and de, where x = x~, represent Cauchy horizons. The singularities 
at hb, ga, df and ce are timelike. The curves ac and bd represent the spatial infinity (x — oo), which is timelike. To have a 
geodesically maximal spacetime, the extension in the vertical directions should be continuous to infinity. 



Case D.2.ii) \q\ = q c : In this case, we find 



/(*) = - *o) V + 2x x + 3x1), (4-46) 



3x 



where xq now is given by 



36 \ 



1/3 



*o=(^ • (4-47) 



Thus, f(x) > holds for all x < 0. As a result, in the region x < the singularity at x = is naked, and the 
corresponding Penrose diagram is that of Fig. 7. 

In the region x > 0, f(x) is also positive, except at the point x — x , where f(x ) = 0. See Fig. 9, and notice that 
now we have x + = x_ = x . This is quite similar to the extreme case of the RN solution [18]. Integrating Eq.(4.7), 
we obtain 

x 1 , (x — x ) 2 

x = -— -r, : r + -tt^ — In ■ 



2\A\x q {x-xq) 3|A|a;o x 2 + 2x Q x + 3x% 
7V2 



x + Xq 



+ — nn arctg + x*, (4.48) 

12|A|a;o \ V2x / 

where Xq is an integration constant. In the following, without loss of generality, we shall choose Xq such that x*(0) = 0. 
On the other hand, from Eqs.(4.6) and (4.46) we find 
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(4.49) 



That is, the surface gravity now is zero. As a result, the usual method of the extension given by Eqs.(4.7)-(4.9) is not 
applicable to this case, and we have to consider other possibilities. Following Carter [19], let us define the two null 
coordinates v and w by 



v = tan 1 (t + x*), w — tan 1 (t — x*), 
where — tt/2 < v < 3tt/2 7 and — 3n/2 < w < tt/2. In terms of v and w, the metric reads 

m 



ds l = — 



-dvdw + x (dy + dz ) 



Using the relation 



sin 2 (w — w) 



4(x*) 2 



(4.50) 



(4.51) 



(4.52) 



one can show that the metric coefficients of Eq.(4.51) become regular across the hypersurface x — xq. Thus, it 
represents an extension. As a matter of fact, one can show that this extension is analytic and the extended spacetime 
is given by regions I and III in Fig. 11, which is quite similar to the extreme RN black hole [18]. The only difference 
is that in the present case the hypersurface x = +oo is time-like, since as x — > +00 we now have x* — > finite. Keeping 
this in mind, we can deduce all the properties of this solution from the extreme RN black hole. For example, the 
horizons x = x represent event horizons. The region x a < x < +00 is mapped to region /, while the one < x < x 
to region III. The t-coordinate is future-directed in region /. The singularity at x = 0, represented by the vertical 
double- lines in Fig. 11, is time- like. From the diagram we can see that the extension should be continuous to infinity 
in the vertical directions, in order to get a geodesically complete spacetime. 




FIG. 11. The Penrose diagram for the solutions given by Eq.(4.46) in the text with \q\ = q c and in the region x > 0. The 
spacetime singularity at x = 0, represented by the vertical line ac, is timelike. The surfaces ae,ab,cb, and cd, where x = xo, 
all represent event horizons. The hypersurfaces x = +00 now are timelike, given by the curves bd and eb. The extension in the 
vertical directions is continuous to infinity. 
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Case D.2.iii) \q\ > q c : In this case it can be shown that f(x) > holds for any x. As a result, the singularity at 
x = is naked, and the corresponding Penrose diagram is given by Fig. 7. 



V. TYPE II FLUID IN D-DIMENSIONAL SPACETIMES 



Now using the gauge freedom of Eq.(2.2), we shall set 

gn(x c )=0, s = r, 



(5.1) 



where x a — {v,r}. These coordinates are usually called Eddington-Finkelstein coordinates, in terms of which the 
metric (2.1) takes the form, 



ds 2 = -e^ v ' r Uv (/(«, r)e^ v ' r Uv + 2edr) + r 2 h AB (x c ) dx A dx B , 



(5.2) 



where e = ±1. When e = 1, the radial coordinate r increases towards the future along a ray v = Const. When 
e = — 1, the radial coordinate r decreases towards the future along a ray v — Const. 
From Eq.(2.8) it can be shown that 



{2) R h 
A(v,r) 
V a V b r : 



[2f^ rr + l rr + (2/Vv + 3/, r ) tA - 2eV,t 
-\fe 2 ^ ^ftr + f,r + ee-^y° a 8° b 

(2/V,r + f,r) (S'JI + S 1 ^) - ijj^jl 



Ur = /</v + /, r . 
Inserting the above expressions into Eq.(2.7), we obtain 



V 



Vtrr + f,rr + (2/^ + 3/, r ) V,r + (2/V>,r + f,r) 



rf 



1 

2 £e 



D - 2 



2/VVr + f,rr + (2/V\r + 3/, r ) 1p, r + 



D-2 



m >r + f >r ) 



- tp,r 



{D) Rab = {k D -r (/W + l r ) — (D — 3)f] h AB . 
On the other hand, introducing two null vectors 1^ and n M via the relations 

Z A ? A = = n A n A , ; A n A = -1, 
we find that the energy-momentum tensor for a Type 77 fluid can be written as [20,21] 

{d) t^ = fxy v + (p + p) (i 

Then, from the Einstein field equations we find that 

{D) R^ = 8nG D j/ii„/„ + (p + P) (!„«„ + J„n M ) + -JfL-g^ . 



(5.3) 



(5.4) 

(5.5) 

(5.6) 
(5.7) 



(5.8) 
(5.9) 

(5.10) 



From Eqs.(5.6) and (5.10) we obtain = 0, that is, ip(v, r) — ijj(v). However, by introducing another null coordinate 
v = J e^^dv, we can always set, without loss of the generality, 
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ip(v, r) = 0, ( (Z3) G rr = 87rG D (D) T rr ) . 

On the other hand, from the rest of the Einstein field equations, we find that 

e(D-2) 



(5.11) 



' f,V ) 



16nG D r 

„■',. 2 ~ (k D -rf, r -(D-3)f), 



P = 



16irG D r 2 
1 

16irG D 



f,rr + ^—^-U — [kn -(D-3) f] } . 



(5.12) 



To see that the energy-momentum tensor given by Eq.(5.9) indeed represents a Type II fluid, let us first introduce 
the following unit vectors, 



i n V 



V2 ' 



(5.13) 



where n a p = diag.{— 1, 1, 1, 1}, and E? A ^ (A, B = 2, 3, D — 2) consist of an orthogonal base in the space of Kab, 
that is, 



J (A)^(B) 

Then, projecting T M „ into this base, we find that 



E ?A) E (B)h-CD — Sab- 



(5.14) 





i(M-2p) 





. 


• °\ 


-\n 





. 


. 








p 


. 


. 











p . 


. 
. 


V o 








. 


. pj 



(5.15) 



which is exactly the Type II fluid, according to the classification given in [18]. Therefore, for any given function 
f(v,r) with ip = 0, the corresponding energy-momentum tensor of the metric (5.2) can be always written in the form 
of a Type II fluid. However, not for any choice of f(v, r) the solution is physically acceptable. In fact, it must satisfy 
some (geometrical and physical) conditions, such as, the energy conditions [18]. Following [21], we write f(v, r) in the 
form 



,, > 2m(v,r) 
f(v,r) = k- rD _ 3 , 

oo 

m(v,r) = a n (v)r n . 



(5.16) 



Then, Eq.(5.12) becomes 



e(D-2) ~ 

TL— — OO 



8irG D r L 



D - 2 °° 

^2 S nfl »W r " _1 ' 



P — 8^^ E »(»-lK(«)r" 



(5.17) 



where a prime denotes the ordinary differentiation with respect to the indicated argument. Clearly, by properly 
choosing the coefficients a n (v), we can get many exact solutions. For example, for the choice, 
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a n (v) 



m(v), 



n = 0, 



-\q 2 {v), n =-(£>- 3), 



2 

±A, 



n = D-l, 
otherwise, 



we find that 



r 2(D-3) 



and the corresponding energy-momentum tensor can be written as 



with 



(D-2)(£>-3) 



1/2 



n-2 W U M U "/ ' 



<D ~ 2) 
87rGi5r B - 2 



m(^)-g(w)-oT3 



(5.18) 



(5.19) 



(5.20) 



(5.21) 



When k = 1 and D = 4 the above solutions reduce to the ones found in [21]. Lately, these solutions were generalized 
to high dimensional spherical spacetimes [22]. 



VI. FORMATION OF TOPOLOGICAL BLACK HOLES FROM GRAVITATIONAL COLLAPSE OF A 

TYPE II FLUID 



To consider gravitational collapse of the Type II fluid found in the last section let us choose e = — 1. Then, metric 
(5.2) with ip = becomes 



ds 2 = -f(v, r)dv 2 + 2dvdr + r 2 h AB (x c ) dx A dx B . 
To the present purpose, let us consider the solutions given by Eqs.(5.19) - (5.21) with 



(6.1) 





r a^- 3 , 


V > Vq, 


m(v) = < 




< v < v n 






v < 0, 




r svo°- 3 , 




q(v) = < 


5v D ~\ 


< v < vq 






v < 0, 



(6.2) 



where A and <5 are arbitrary real constants, subject to A > 0. Clearly, the solutions represent a charged null dust fluid 
moving on a de Sitter or anti-de Sitter background in the region < v < vq, depending on the sign of A. When A = 
the spacetime is self-similar, while when A ^ it is only asymptotically self-similar, (v, r) — ► (0, 0). The corresponding 
energy-momentum tensor is given by Eq.(5.20) with 



(D) F 



( 



-2)(D-3) \ 1/2 6vo°-3 



M = 



(D-2)(g-3) 
2 J 
0, 

(D-2)(D-3)8 2 y L 
8irG D r 2 

0, 



1/2 



r D-2 
.,D-3 



- <W . < v < vo, 
v < 0, 



D-3 



,D-3 



) , < V < Vq, 

otherwise, 



where 



V = 



Vc = 



S 2 



l/(D-3) 



(6.3) 



(6.4) 
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From Eqs.(5.4)-(5.7) we can see that the Ricci tensor contains only the first order derivatives of / with respect to v. 
Thus, there would be no matter shell to appear on a hypersurface v = Const., as long as f(v,r) is continuous across 
this surface. Clearly this is the case for the choice of Eq.(6.2) crossing the hypersurfaces v = and v = vq. 

When k = 1, A = = 5, and X D ~ 2 = S D ~ 2 , the corresponding solution is the Vaidya solution in D-dimensional 
spherically symmetric spacetime [22], and was studied independently in [23,24]. In particular, da Rocha showed that 
when A > A c the collapse always forms black holes, and when A < A c the collapse always forms naked singularities, 
where 



(D-3) 



D-3 



[2{D-2)\ D - 



(6.5) 



This result can be easily generalized to the case where A ^ 0, 5 — and X D ~ 2 has a topology rather than S D ~ 2 . As 
a matter of fact, the case with k = 1, A>0, S = and X D ~ 2 = S D ~ 2 was already studied in [25], from which we 
can see that Eq.(6.4) is also valid for any A and other topologies for X D ~ 2 but still with k = 1, (5 = 0. 

When k = 1, A = 0, (5^0 and D = 4, Lake and Zannias found that the collapse can form either black holes or 
naked singularities [26]. Ghosh generalized these results to the case where A > [27]. From the analysis given in 
these two articles one can see that the cosmological constant actually has no effects on the final state of the collapse. 
Therefore, the results obtained by Lake and Zannias are actually valid for any A. In addtion, following their analysis 
it is not difficult to be convinced that gravitational collapse of a charged Type II fluid in a higher dimensional de 
Sitter (A > 0) or anti-de Sitter (A < 0) spacetime with k = 1 can form either black holes or naked singularities, 
depending on the choice of A and S. 

Therefore, in the following we shall consider only the cases where k = 0, —1. Let us first consider the case 8 = 0, 
i.e., the collapse of a neutral null dust fluid. In this case from Eq.(5.19) we can see that to have f(v, r) be positive at 
least in some regions of the spacetime, we have to assume that A < 0. Then, Eq.(5.19) can be written as 



f(v,r)=2X 



B [r 2 r 2 ) yDS 

1/2 



(6.6) 



while Eq.(A.17) yields, 



0i6 n = 



2XFG 

„D-1 



D-3 



IAI 




2A 

v < v AH (r), 
v = v A u(r), 
v > v AH (r), 



where 



vah{t) = r j j^j (r 2 - r 2 )| 



l/(D-3) 



(6.7) 



(6.8) 



Thus, in the present case the {D — 2)-surfaces of constant v and r are trapped in the region v > VAuir)- The 
hypersurface v — vah (r) represents an apparent horizon [cf. Fig. 12]. Since the spacetime singularity at r = 
starts to be formed only at the moment v = 0, from Fig. 12 we can see that this singularity is always covered by the 
apparent horizon. This can be seen further by studying "outgoing" null geodesies, which are given by 



dv 



dy |A|r 2 -|fc|-2Ay D - 3 ' 



(6.9) 



At the moment v = if there exist out- going null geodesies from the point (v,r) = (0,0), we can see that the 
singularity formed at that moment will be at least locally naked. The existence of such null geodesies is characterized 
by the existence of positive roots of the equation [28] , 



2Xy D - 2 + \k\y + 2 = 0, 



(6.10) 



where 
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v dv 
yo = lim — = lim — . 

v.r— >o r «,i — »o ar 



(6.11) 



Since all the coefficients of Eq.(6.10) are positive, no positive roots exist. As a result, the collapse will always form 
black holes with non-trivial topology. This generalizes the results obtained in four-dimensional case [10] to the one 
with any dimensions. 




FIG. 12. Collapsing Type II fluid described by Eq.(6.2) in the text with k 
represent an apparent horizon, and the collapse always forms black holes. 



< 0, <5 = 0, and A < 0. The curved line AH 



When <5 ^ 0, Eq.(6.10) should be replaced by 

S 2 yo 2D - 5 - 2\y D - 2 - \k\y -2 = 0. (6.12) 

It is important to note that in the present case the charged fluid satisfies the weak energy condition only in the region 
V < Vci a s we can see from Eq.(6.3). On the hypersurface y — y c the energy density becomes zero. Afterwards, the 
Lorentz force will push the fluid particles to move outwards [29]. As a result, the particles actually cannot enter into 
the region y > y c [cf. Fig. 13]. Now to see if the spacetime singularity formed at (v,r) = (0,0) is naked or not, one 
needs not only to show that a positive root of Eq.(6.12) exists, but also to show that the outgoing null geodesies fall 
inside the region y < y c , where the solution is actually valid. 



II 



^ ^ Singularity 


/y 


>> 


- (t 


) 




/ v 

/ y = y c 

- (a) 
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<y 


c 







FIG. 13. Null geodesies for collapsing Type II fluid described by Eq.(6.2) with 8^0. The dashed line (a) represents an 
outgoing null geodesic with j/q < Vc while the one (6) an outgoing null geodesic with > Vc where denotes the positive 
root of Eq.(6.12) and y c is defined by Eq.(6.4). 
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From Fig. 13 we can see that if an outgoing null geodesic falls into the region y < y c we must have y$ < y c , and if 
an outgoing null geodesic falls into the region y > y c , we must have t/g > y Cl where j/g is the smallest positive root 
of Eq.(6.12). Since in the region y > y c we have fi < and in a realistic model this region should be replaced by 
an out-going charged dust fluid. Thus, in the latter case the singularity at (v,r) = (0,0) should not be considered 
as naked, but in the former case it is. Therefore, to see if the singularity is naked or not now reduces to find out if 
Vo < Vc ° r Vo > Vc- To this end, let us first consider the case k = 0, and define the function G (yo) by 

G(y ) = d 2 y 2D - 5 -2\y D - 2 -2, (k = 0). (6.13) 

From this expression we can see that G" (yo) = has two roots, y^ , given by 

and 

G(vo) =-2, 

_ 4A 2 QD-2)QD-3) 1/(J - S) 
^{Vo) - 5 2 {2D _ b) 2 Vo 2< 2, 

G ^=-p{p) (6 ' 15) 

Then, the curve of G (yo) vs yo must be that given by the line (a) in Fig. 14, from which we can see that 

J/o > Vc (6-16) 

That is, in the present case there does not exist outgoing null geodesies in the region y < y c . So, the singularity 
formed at the point (v, r) = (0, 0) is not naked. 
When k = — 1, setting 

F(y )^S 2 y 2D - 5 -^yo D - 2 -Vo-2, (k = -1), (6.17) 

we find that 



G(y )-F(yo) = y , (6.18) 

and the curve of F (yo) vs yo must be given by the dashed line (b) in Fig. 14. From there we can see clearly that 
Eq.(6.16) also holds for k = —I. 





G(y ) 








(a) 







■/ / (b) 






-\ / / y ° ' y ° 



FIG. 14. (a) The function of G (yo) defined by Eq.(6.13) for k = 0. (b) The function of F (yo) defined by Eq.(6.17) for 
k = -1. 
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Therefore it is concluded that the gravitational collapse of a charged null dust fluid given by Egs.(5.19), (6.1) and 
(6.2) with fc = or k ~ — 1 always forms black holes. This is in contrast to the results obtained in [27]. The reason 
is that in [27] the author didn't consider the energy condition of the charged null dust fluid. 

VII. CONCLUSIONS 

In this paper, we have studied topological black holes and their formation from gravitational collapse of a Type 
II fluid in D-dimensional spacetimes described by the metric Eq.(2.1). After presenting a general R 2 x X D ~ 2 
decomposition in Sec. 77, we have re-derived the charged solutions coupled with a cosmological constant, given by 
Eq.(3.21), in Sec. III. From their derivations we can see that these solutions are the most general ones for the 
spacetimes described by Eq.(2.1). 

In Sec. IV, we have systematically studied the global structure of the spacetime for kr> = and D = 4, that is, four- 
dimensional spacetimes with flat topology of the X 2 sector, which can be a two-dimensional plane, cylinder, Mobius 
band, torus, or Klein bottle, depending on how to identify the two coordinates in X 2 [13]. All the corresponding 
Penrose diagrams have been given. In particular, it has been found that the solution with A < and q = has a black 
hole structure quite similar to the Schwarzschild black hole [cf. Fig. 9]. In the case where A < 0, bq^{), the global 
structure of the corresponding spacetime is quite similar to the Reissner-Nordstrom solution, including the extreme 
case [cf. Figs. 10 and 11]. 

In Sec. V, all the solutions of a Type II fluid have been found, while in Sec. VI the gravitational collapse of such 
a fluid has been studied. When k = 1 the collapse in a de Sitter or anti-de Sitter background can form either black 
holes or naked singularities, but when k — or k = — lit always form black holes. Therefore, all the black hole 
solutions with different topologies found in Sec. IV can be realized from the gravitational collapse of a type II fluid 
in D-dimensional spacetimes. 
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APPENDIX: TRAPPED SURFACES AND APPARENT HORIZONS 

The concept of trapped surface was originally from Penrose [30] , who defined it as a compact spatial two-surface S 
in a four-dimensional spacetime, on which 9+9-\ s > 0, where 9± denote the expansions in the future-pointing null 
directions orthogonal to S, and the spacetime is assumed to be time-orientable, so that "future" and "past" can be 
assigned consistently. One may then define a past trapped surface by 9±\ s > 0, and a future trapped surface by 
9 ± \ s <0. 

Recently, Hayward [31] generalized the above definition to the four-dimensional cylindrical spacetimes where the 
two-surface S is not compact but an infinitely long cylinder, and call it trapped, marginal or untrapped, according to 
where is timelike, null or spacelike, where a is defined by 

a = (\d z -d z \-\d v -d v \) 1/2 , (A.i) 

with £( z ) = d z and = d v being the two killing vector in the cylindrical spacetimes. Lately, one of the present 
authors showed that Hayward's definition is consistent with the one defined by the expansions of the two null directions 
orthogonal to the two-cylinder [32]. 

In this Appendix, we shall generalize the above definition to the D-dimensional case described by the metric (5.2). 
Now let us consider the (D-2)-surface, S, of constant v and r. When k = 1 it could be a (D — 2)-sphere, although 
other topologies can also exist. When k ^ 1 it can be compact or non-compact, depending on how to identify the 
coordinates in S. For more details, we refer readers to [13,33] and references therein. 

To calculate the expansions of the two null directions orthogonal to the (D — 2)-surface S of constant x a in the 
metric (5.2), it is found convenient to introduce two null coordinates v and u via the relations 

du = G{v, r) (fe^dv + 2edr) , 

dv = F{v)dv, (A.2) 
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or inversely 

dv = ^dv, (A.3) 
F 

where G(v,r) is determined by the intcgrability condition u yVr = u, rv , and F(v) is an arbitrary function of v only. 
Without loss of generality we shall assume that they are all strictly positive, 

F(v)>0, G(v,r)>0. (A.4) 

Then, in terms of u and v the metric (5.2) takes the form 

ds 2 = -2e 2rT dudu + r 2 h AB (x c ) dx A dx B , (A.5) 

where a and r are now the functions of u and v via Eq.(A.2) and 

a=^-ln(2FG)}. (A.6) 

Clearly, the metric (A.5) is invariant under the coordinate transformations 

u = u(u) 7 v — v(v). (A-7) 

Using this gauge freedom, we shall assume that metric (A.5) is free of coordinate singularities. Then, it can be shown 
that 

WR ab = -2* m (sy b + SX) , 

V a V 6 r = r, ab - 2 {a^r^Jl + a^r^jl) . (A.8) 
Substituting Eq.(A.8) into Eq.(2.7) we find that 

(D) d _ D ~ 2 (r 0(7 -r -) 

1L uu r ' ,UU - u - u ' ' 

(D) R _ o a D ~ 2 r 
r 

(D) R __ _ ( r 2rr -r -) 

r 

{d) Rab = {k D + 2e- 2a [rr iU v + (D - 3)r iU r,v}} h AB . (A.9) 
On the other hand, introducing two null vectors and by 

one can see that the two null vectors are future directed, and orthogonal to the (D — 2)-surface, S, of constant u and 
v (or constant of v and r). In addition, each of them defines an affinely parametrized null geodesic congruence, since 
now the following holds 

■^-n" = n^, v n v = 0, (A.ll) 
L>o 

where a semicolon ";" denotes the covariant derivative with respect to g^, and A and <5 the affine parameters along 
the null rays defined, respectively, by and (It should be noted that the symbol 8 used here should not be 
confused with the one used in Eq.(6.2) for the charge density.). In particular, defines the one moving along the 
null hypcrsurfaccs u = Const., while defines the one moving along the null hypersurfaces v — Const. Then, the 
expansions of these null geodesies are defined by [18], 
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Thus, we have 



~9^n a; , = e-^ r f. (A.12) 



OlOn = e-^-^f = -^e- 2ff W*. (A.13) 



It should be noted that the affine parameter A (or 5) is unique only up to a function f (u) (or g _1 (u)), which 
is a constant along each curve u = Const, (or v — Const.) [18]. In fact, A = A/ f(u) (5 — S/g(v)) is another affine 
parameter and the corresponding tangent vectors are 

r„ = /(«)** n^=g(v)Sl, (A. 14) 

and the corresponding expansions are given by 

d T =f{u)8 u 6 n = g(v)6 n . (A.15) 

However, since along each curve u = Const, (or v = Const.) the function f(u) (or g(v)) is constant, this does not 
affect our definition of trapped surfaces in terms of the expansions. Thus, without loss of generality, in the following 
we consider only the expressions given by Eq.(A.12). 

Once we have the expansions, following Penrose we can define that a (D — 2)-surface, S, of constant u and v is 
trapped, marginally trapped, or untrapped, according to whether 6i9 n > 0, 9iO n — Q, or9iO n < 0. An apparent horizon, 
or trapping horizon in Hayward's terminology [31], is defined as a hypersurface foliated by marginally trapped surfaces. 

Since e~ 2a is regular, except at some points or surfaces on which the spacetime is singular, from Eq.(A.13) we can 
see that trapped, marginally trapped, or untrapped surfaces can be also defined according to whether r Q is timelike, 
null, or spacelike. 

On the other hand, from Eq.(A.3) we find that 

* = _L — = -lL </> 
du " 2G" dv ~ 2F 6 ■ 

Inserting the above expressions into Eq.(A.12) we obtain 

a — ^ r, a Pl — ~' 2 <r r > ! " 

D r 



(A.16) 



jf_ 

2rF 

D r 

e „-2<r 



2rG 



FC 

~fe~^. (A.17) 
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